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EVALUATION OF MASSLESS-SPRING MODELING OF 
SUSPENSION- LINE ELASTICITY DURING THE 
PARACHUTE UNFURLING PROCESS 

By Lamont R. Poole and Earle K. Huckins III 
Langley Research Center 

SUMMARY 

Mathematical modeling of elasticity in parachute suspension lines during the 
unfurling process was investigated. The massless- spring approximation was evaluated 
in depth. Equations which govern the elastic response of the suspension lines for this 
basic model were combined with equations governing the motion of the vehicle and the 
deployment bag. This set of governing equations was integrated numerically and the 
results were compared with flight test data. Agreement was satisfactory except during 
periods of rapid fluctuations in line tension. Lack of agreement during these periods 
was attributed primarily to omission of wave mechanics in the suspension-line elasticity 
model under evaluation. 


INTRODUCTION 

Parachute deployment is defined as the process of extending the parachute canopy 
and suspension lines to a point at which satisfactory inflation can occur. Interest in 
decelerator systems for planetary entry vehicles has stimulated research in deployment 
dynamics which has resulted in significant improvements over earlier empirical methods 
(refs. 1 and 2). A primary need for simulating the lines- first deployment process is to 
determine the motion of the deployment bag relative to the towing vehicle. Knowledge of 
sequence of events and deployment times is important from a mission-analysis stand- 
point. A greater understanding of the process by which the packed parachute is unfurled 
is important in deployment system design. An accurate prediction of the tension 
developed in the suspension lines during lines-first deployment is important from a 
standpoint of vehicle motions. 

An analysis of the lines-first deployment process for an inelastic parachute was 
presented in reference 3. The analysis presented in reference 4 indicated that the 
effects of suspension- line elasticity can be significant during deployment. The effects of 
suspension- line elasticity during the parachute inflation process have been studied pre- 
viously (ref. 5) for a simple two-body spring-mass system. Mathematical modeling of 



elastic suspension lines during the unfurling process is discussed in appendix A. The 
purpose of the present paper is to evaluate the use of a simple massless -spring model in 
the analysis of lines-first deployment dynamics. 

In the present paper, equations are developed which, for a lines-first forced- 
ejection deployment, govern the planar motion of the vehicle, the linear motion of the 
deployment bag relative to the vehicle, and the unfurling process. In the development, 
the parachute suspension lines are represented by a variable-length massless spring. 

The usefulness of this mathematical model is evaluated by employing it in the calculation 
of deployment motions and loads of two disk-gap-band parachute deployment flight tests. 
Uncertainty in the mechanisms by which the packed parachute is unfurled from the deploy- 
ment bag is accounted for by using both precise and averaged linear-mass-density distri- 
butions for each of the parachutes. 


SYMBOLS 


A reference area, meters 

c wave propagation velocity of suspension lines, meters/second 

C damping coefficient of a single suspension line, newton- seconds 

Cn drag coefficient, 

u q A 

oo 

Fj-j t aerodynamic drag per unit length of suspension lines, newtons/meter 

F re unfurling resistance force, newtons 

g acceleration due to gravity, meters/second^ 

h altitude, meters 

K' specific secant modulus of a single suspension line, newtons 

sec 

l general spatial coordinate along unstressed suspension lines, measured from 

vehicle attachment point, meters 

unfurled length (or, equivalently, spatial coordinate of parachute mass ele- 
ment exiting mouth of deployment bag), meters 
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total unstressed length of extended parachute, meters 


total unstressed length of suspension lines, meters 
mass, kilograms 

linear mass density of general element (parachute) exiting mouth of 
deployment bag, kilograms/meter 

linear mass density of a single suspension line, kilograms/meter 

number of suspension lines 

free-stream dynamic pressure, newtons/meter^ 
time, seconds 
tension, newtons 

unfurling velocity or rate at which parachute exits mouth of deployment bag, 
meters/second 

velocity, meters/second 

ejection velocity, meters/second 

deployed distance relative to vehicle, meters 

vehicle flight-path angle, degrees 

displacement of any cross section of suspension lines from its unstressed 
position, meters 

strain, meters/meter 

wake parameter (ratio of drag coefficient of deployment bag in vehicle wake 
to that in free stream) 


P to free-stream density, kilograms/meter^ 

co effective loading frequency, radians/second 

Subscripts: 

b deployment bag 

B deployment-bag mouth 

dyn dynamic 

e deployment bag plus its instantaneous contents 

o conditions at beginning of deployment 

s static 

v vehicle 

V suspension-line attachment point at vehicle 

Dots over symbols denote time derivatives. 

ANALYSIS 

In the analysis to follow, equations which govern the dynamics of lines-first deploy- 
ment of a parachute having elastic suspension lines of the massless-spring class are 
developed. Equations which govern the motion of the vehicle and the deployment bag are 
presented. An expression for the rate at which the parachute exits the mouth of the 
deployment bag (the unfurling rate) is given. In addition, equations are presented by 
which the elastic behavior of the suspension lines can be approximated. 

A typical vehicle-parachute configuration during a deployment of the lines-first 
type (ref. 4) is shown in figure 1. A sketch of the forces affecting the dynamics of the 
deployment process is given in figure 2. For the purposes of the present analysis, the 
surface of the earth is assumed flat, and surface- relative accelerations are considered 
inertial. 
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Motion of the Vehicle 


The motion of the vehicle is influenced by aerodynamic forces, tension in the para- 
chute suspension lines, and gravitational attraction. Assume that 

(1) The vehicle can be represented by a mass particle (vehicle moments of inertia 
are neglected). 

(2) The suspension- line tension at the vehicle is parallel to the vehicle relative 

wind. 

(3) The vehicle is nonlifting. 

(4) Winds are negligible. 

Under these assumptions, the longitudinal acceleration of the vehicle is given by 


v v = - 


( CdA K 


+ T, 


m. 


+ g sin y 


( 1 ) 


where 


inertial acceleration 


C D 


drag coefficient 


A 

q 


reference area 


1 2 

free- stream dynamic pressure Iq^ = ^p^v v 


free- stream density 


tension at suspension- line attachment point on vehicle 


m 


mass 


g 


acceleration due to gravity 


vehicle flight-path angle 


and subscript 


vehicle 
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In order to completely specify the earth- relative planar motion of the vehicle, the 
following additional trajectory equations are required: 


v v sin y 

(2) 

g cos y 
v v 

(3) 


where h v is the time rate of change of the altitude of the vehicle. 


Motion of the Deployment Bag 


The motion of the deployment bag can be closely approximated by introducing the 
following additional assumptions: 

(1) The deployment bag and its instantaneous contents can be represented by a 
variable-mass particle. (Moments of inertia are neglected.) 

(2) The velocity vector of the deployment bag relative to the vehicle is parallel to 
the relative-wind velocity vector of the vehicle. 


(3) The reduction in dynamic pressure acting on the deployment bag due to the 
motion of the bag relative to the vehicle is neglected. 


Under these assumptions, the longitudinal acceleration of the deployment bag relative to 
the inertial frame is given by 


v b - 


' - F re 


m £ 


+ g sin y 


(4) 


where 

r\ wake parameter 

F re unfurling resistance force 

and subscripts 
b deployment bag 

e deployment bag plus its instantaneous contents 


6 



In addition, the velocity of the deployment bag relative to the vehicle can be written 
as follows: 


Xu = V„ - V, 


(5) 


Unfurling Rate 

A relationship between the tension in the unfurling parachute at the mouth of the 
deployment bag and the unfurling rate was given in reference 4, and can be written as 

Tg = m'u^ + F re (6) 


where 


Tg tension in unfurling parachute at mouth of deployment bag 

m T linear mass density of parachute exiting mouth of deployment bag 

u rate at which parachute exits mouth of deployment bag (unfurling rate) 

and 

« 

u = Z B 

or, by using equation (6), 



(?) 


(8) 


Suspension- Line Elasticity 

In order to completely define the motion of the deployment process, it is necessary 
to formulate a relationship between the suspension-line tension at the vehicle, the tension 
at the mouth of the deployment bag, the motion of the deployment bag relative to the vehi- 
cle, the unfurling rate, and the unstressed length of the unfurled parachute. This relation- 
ship is furnished by a mathematical model of the suspension- line elasticity. 

Mathematical modeling of elastic suspension lines is discussed in appendix A. It 
is shown that, under the following assumptions, the elastic suspension lines can be approx- 
imated by a massless spring of variable length: 

(1) Throughout the deployment process, the suspension lines are under tension. 
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(2) Aerodynamic drag on the unfurled suspension lines is negligible. 

(3) Suspension- line tension is independent of strain rate. (Viscous damping is 
neglected.) 

(4) Linear mass density of the suspension lines is negligible. 

(5) The effective loading frequency is much less than the effective propagation 
frequency. (See appendix A.) 

The tension throughout the suspension lines is constant at any particular time: 

T = Ty = T b (9) 

where T is the tension at any arbitrary spatial location along the unfurled suspension 
lines and is given by 

T.n sl K’ ec e (10) 

where 

n s i number of suspension lines 

K' specific secant modulus of a single suspension line 

e strain in suspension lines 

The suspension- line strain is by definition 


X B ' 

l B 

L sl) 

(11) 

X B " ^B 
L sl 

(Lsl < l B < L p) 

I (12) 


where 

Xg deployed distance (distance from suspension- line attachment point on vehicle 

to deployment-bag mouth) (Note: Xg = x^ 

L s i total unstressed length of suspension lines 

Lp total unstressed length of extended parachute 
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RESULTS AND DISCUSSION 


A mathematical model to describe precisely the dynamics of the entire parachute 
deployment sequence would necessarily involve direct solution of the partial differential 
equations governing the wave motion of the suspension lines. In the present paper, the 
consequences of neglecting wave mechanics were studied by evaluating a massless-spring 
model of suspension- line elasticity. 

In order to evaluate the accuracy obtained by using this simple model, it was used 
to calculate the deployment motions for two flight tests of disk-gap-band parachutes: the 
second balloon- launched flight test of the NASA Planetary Entry Parachute Program 
(ref. 6) and the first flight test of the NASA Supersonic Planetary Entry Decelerator Pro- 
gram (ref. 7). The governing equations, in conjunction with the auxiliary expressions 
described in appendix B, were solved numerically on a digital computer, using a fifth- 
order Runge-Kutta integration technique, for appropriate values of parameters and coef- 
ficients. In order to allow for uncertainties in the mechanisms by which the packed 
parachute is unfurled from the deployment bag, two linear mass distributions for each 
parachute were employed in the solution: a precise mass distribution (of the form shown 
in sketch B1 of appendix B) and a distribution consisting of averaged values of linear 
mass density (average of the precise values over a characteristic parachute folding 
length, the deployment-bag diameter). 

Calculated histories of unfurled length for the two flight tests are presented in 
figure 3 and compared with test data points, for which occurrence times have been 
adjusted to account for wave propagation times. As neither of the calculated histories 
was affected significantly by the alternate mass distribution, only one such history is pre- 
sented for each flight. The flight test data points presented correspond to mortar fire, 
line stretch, and estimated bag strip (ref. 8). The calculated histories show good agree- 
ment with flight test data. 

Calculated histories of tension for the two flight tests are presented and compared 
with test data in figure 4. Again, neither of the calculated histories was significantly 
affected by the alternate mass distribution, and only one tension history is presented for 
each flight. The initial large peak in each flight test data history corresponds to the 
impact of a tensiometer, which was included in the flight system in order to measure 
tension levels. This occurrence was not considered in the present mathematical model. 
Agreement is acceptable except during the periods corresponding to the unfurling of the 
band portions of the respective canopies. Lack of agreement during these periods of 
rapid load fluctuations is attributed to omission of wave mechanics in the suspension- 
line elasticity model. 
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Calculated histories of the deployment rate and unfurling rate using both precise 
and averaged linear mass distributions are presented in figure 5 for the flight test of ref- 
erence 6 and in figure 6 for the flight test of reference 7. Similar trends are exhibited in 
each figure. The deployment-bag motion is quite smooth in nature, independent of the 
choice of linear mass distribution. In general, the unfurling rate history closely follows 
the deployment rate history. However, sharp decreases in the unfurling rate occur when 
sudden increases in the linear mass density of the unfurling parachute are encountered; 
conversely, sharp increases in unfurling rate occur when sudden decreases in linear 
mass density are encountered. The effect of averaging the linear mass density over a 
characteristic length is virtually to eliminate the narrow "spikes" in the precise mass 
distribution which are representative of cloth overlapping and tape reinforcement in the 
parachute canopy. As a consequence, the magnitudes of the decreases in unfurling rate 
are reduced when the averaged mass distribution is used. The increases in unfurling 
rate, which occur as the gap and vent portions of the canopy unfurl from the bag, are not 
affected significantly by the alternate mass distribution. As the gap emerges from the 
bag, the unfurling rate actually exceeds the ejection velocity of the respective flight tests. 
In consideration of the relatively poor agreement of calculated tension histories with 
flight test data and the absence of data concerning unfurling rate, some uncertainty exists 
as to the magnitudes of the calculated maxima and minima in unfurling rate. However, 
the trends exhibited in the calculated histories are representative of the physical 
situation. 


CONCLUSIONS 

Massless-spring modeling of suspension- line elasticity during the lines-first para- 
chute unfurling process has been studied. Based on numerical solution of the governing 
equations for two disk-gap-band parachute deployment flight tests, the following conclu- 
sions can be made in evaluation of the present model: 

1. Calculated histories of unfurled length were not affected significantly by either 
the elastic model or the alternate parachute linear mass distributions and agreed very 
well with flight test data. 

2. Calculated tension histories agreed poorly with flight test data during periods 
corresponding to the unfurling of the band portion of the canopies. Acceptable agreement 
was obtained during periods corresponding to suspension- line and disk unfurling. Ten- 
sion histories were not affected significantly by the alternate linear mass distributions. 

3. Calculated histories of unfurling rate exhibited fluctuation about the smooth 
deployment rate histories, generally as a function of the selected linear mass distribu- 
tion. Inasmuch as the local maxima in the unfurling rate histories were not affected 
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significantly by the alternate mass distributions, there exists the possibility that 
unfurling rates of greater magnitude than the flight- test ejection velocities were experi- 
enced during the periods corresponding to unfurling of the gap portion of the respective 
canopies. 

4. In light of the high-frequency loading experienced during canopy unfurling, a 
formulation which considers the wave mechanics of the suspension lines will probably 
be required in order to reduce the remaining inaccuracies and to clear up existing 
uncertainties in the unfurling process. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 3, 1972. 
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APPENDIX A 


MATHEMATICAL MODELING OF PARACHUTE 
SUSPENSION LINES 

Mathematical models of parachute suspension lines can vary greatly in complexity. 
The motion of the suspension lines during deployment can be determined in detail by 
writing Newton's second law for an element of the lines, formulating an elastic model 
for the element, and integrating the resulting set of coupled nonlinear partial differential 
equations. The results obtained would completely describe the lateral and longitudinal 
motions of the suspension lines. 

However, this degree of complexity is not usually required since, in many cases, 
the deployment motion is approximately one-dimensional. Under the assumption of linear 
motion, the governing equations reduce to a single nonlinear partial differential equation 
which describes the longitudinal wave motion in the suspension lines. This substantial 
reduction in the complexity of the analysis introduces relatively small inaccuracies in 
the results. Linear motion will be assumed in the analysis to follow. 

The elastic state of the suspension lines can be represented by use of the coordi- 
nates and variables shown in sketch Al. 



Vehicle 


Sketch Al 
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APPENDIX A - Continued 


The longitudinal forces acting on an element of the suspension lines are shown as 
follows: 


where 


T -t- 


dZ 


T + (fz + n si m sis sin y + f d') dl 


tension in element 


m sl 


linear mass density of a single suspension line 


n sl 


number of suspension lines 


acceleration due to gravity 


vehicle flight-path angle 


'D 


aerodynamic drag per unit length of suspension lines 


Assuming linear motion, Newton's second law for the longitudinal motion of the line 
element can be written as 


8T 


~qT + n sl m slS sin y + F d' = n sl m sl “2 - v v 


. 2 6 


ar 


(Al) 


where v y is the inertial acceleration of the vehicle given by 


where 


Vv = - 


( C D A ) v^o 


■V 


m 


v 


+ g sm 7 


(A2) 


'D 


drag coefficient 


reference area 


free-stream dynamic pressure 
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APPENDIX A - Continued 


and subscript 


vehicle 


The governing equation for the motion of the element can be obtained by substituting 
equation (A2) in equation (Al). The result can be written 


9 2 6 


9T 


f d' 


( c D A b. 


+ T V 

OO V 


8t 


2 Hgirngi dl n sl m sl 


m. 


(A3) 


In order to define the motion completely, it is necessary to determine a func- 
tional relationship between the tension, the displacement 6, and its derivatives. The 
tension in the line element is a result of static tension T s and tension due to dynamic 
effects TdynJ that is, 


T T s + Ttfyn 


(A4) 


The static tension can be related to the strain in the following manner: 


T s = n S l K sec e 


(A5) 


where 


, . . 96 

strain ( e = 


K sec 


specific secant modulus of a single suspension line 


n sl 


number of suspension lines 


Note that, in general, the secant modulus is a function of the strain and the strain history; 
and, as such, it can include the effects of nonlinear elasticity, mechanical hysteresis, and 
permanent set. 

In a similar manner, the tension due to dynamic effects can be related to the strain 
rate as follows: 


L dyn 


= n sl c < 


(A6) 
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APPENDIX A - Continued 


where C is the damping coefficient of a single suspension line. In general, the damping 
coefficient is a function of both the strain and the strain rate. This nonlinear function of 
strain rate accounts for elastic hysteresis (the elastic aftereffect) and viscous (or rate) 
damping. 


The tension gradient in equation (A3) can now be written in terms of the displace- 
ment gradients and their time derivatives by use of equations (A5) and (A6). The result 
is 


9T 

U ~ n sl 


K' 3 2 5 , 9K sec 
sec *1 



(A7) 


Initial conditions on the motion are as follows : 

(1) The initial displacement profile of the suspension lines. 

(2) The initial displacement rate profile of the suspension lines. 

In addition, the dynamic response is subject to the following boundary conditions: 

(1) The displacement at the vehicle 

6| =0 (A 8) 

7=0 

(2) The displacement rate at the mouth of the deployment bag 


as 

at 


l=lB 


x b - u 


where 


spatial coordinate of mass element exiting deployment-bag mouth 
unfurling velocity 
velocity of bag relative to vehicle 

(3) The tension at the mouth of the deployment bag 
T B = m'u 2 + F re 


u 

*b 


(A9) 


(A10) 
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APPENDIX A - Continued 


where 


Tg tension in element unfurling from deployment bag 

m' linear mass density of unfurling element 

F re unfurling resistance 


A simultaneous solution to equations (A3) and (A7) subject to the initial and boundary con- 
ditions listed above will completely define the elastic state of the suspension lines. 

In most cases, a solution of the complete system of equations can only be obtained 
numerically. A suggested procedure is to approximate partial derivatives by difference 
quotients formed from the finite difference in elastic states of adjacent elements. Finite 
difference schemes for solution of partial differential equations are discussed in many 
textbooks on numerical analysis. (For example, see ref. 9.) 

The governing equations can be simplified considerably by neglecting damping and 
assuming that the elastic and mass properties are the same at all cross sections along 
the suspension lines. Under these assumptions, equations (A3) and (A7) can be combined 
and reduced to 

a2e K sec 9^5 f d' (An) 

9t 2 m sl a l 2 n sl m sl m v 


For linearly elastic suspension lines (Kg ec = Constant), the wave velocity can be 
introduced: 


K T 

„2 _ sec 
m si 


(A 12) 


where c is the longitudinal wave velocity. 

For linearly elastic suspension lines, negligible drag on the unfurling suspension 
line, infinite vehicle mass, and constant m* p equation (All) reduces to the one- 
dimensional wave equation which has a classical solution. In reference 4, the classical 
solution, for a particular set of initial and boundary conditions, was used to predict the 
"snatch force" which occurs as the canopy skirt emerges from the deployment bag. 

Variations in the linear mass density of the material exiting the mouth of the deploy- 
ment bag result in load fluctuations in the suspension lines. During events such as line 
stretch, the effective frequency of these load fluctuations can be relatively high. As a 
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APPENDIX A - Continued 

result, a loading cycle may be completed before the response wave has traversed the 
entire length of the suspension lines. In order to accurately model the suspension 
lines during such an event, wave mechanics must be considered. However, during 
periods in which the effective loading frequency is much less than the effective propaga- 
tion frequency, that is, 

w « (A13) 

L sl 


where 


co effective loading frequency 

(2 

c wave velocity I c = 

L s i length of unfurled suspension lines 


K sec\ 

m sl/ 


the wave velocity may be considered infinite. Infinite wave velocity tmplies either 
inextensible suspension lines ^Kg ec — or massless suspension lines (ni' s i — 0 ). 

For inextensible suspension lines, the displacement of the lines from the unstressed 
position is identically zero and equation (A3) reduces to 


IT = n si m sl 




+ T V 


m T 


f d 


(A14) 


which can be integrated to give the tension profile in the inextensible lines: 


F F B n sl m sl 



(A 15) 


The inextensible model, simplified further by neglecting aerodynamic drag on the lines, 
was used to study the dynamics of the unfurling process in reference 3. 

For massless suspension lines (note however, that the lines are assumed massless 
with respect to elastic modeling, not with respect to tension generation at the mouth of the 
deployment bag), equation (All) reduces to 


9 2 5 

dl 2 


V 


n sl K sec 


(A 16) 
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APPENDIX A - Concluded 


which can be written 

f£ = - F D’ (A17) 

Assuming that the drag of the unfurled suspension lines is much less than the suspension- 
line tension, equation (A17) can be closely approximated by 

^ = 0 (A 18) 

U 

Therefore, under the assumption of vanishing linear mass density, tension is constant 
throughout the lines at any particular time. 

The dynamic response of suspension lines during the unfurling process can be most 
accurately determined by using a mathematical model which includes wave mechanics. 
However, a numerical solution to the governing partial differential equation is consider- 
ably more complex than a numerical solution to the ordinary differential equations char- 
acteristic of the simpler models. Additional complexities may be introduced in a model 
which includes wave mechanics by the dispersion of the wave velocities between individ- 
ual suspension lines. The purpose of the present analysis is to evaluate a massless- 
spring model of the suspension lines, since this represents the simplest model which 
includes the effects of elasticity. 
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APPENDIX B 


AUXILIARY COMPUTATIONAL EXPRESSIONS 

Based on the assumptions made, the governing equations which were developed in 
the analysis section of this report define the dynamics of the deployment process. How- 
ever, in order to successfully adapt this system of equations to standard numerical inte- 
gration techniques, several additional expressions must be developed. These include 
expressions for the unfurling acceleration and for initial values of the unfurling rate 
and unfurling acceleration. In addition, a technique must be developed to treat disconti- 
nuities in the linear mass density. Several procedures which have been found useful are 
described in this appendix. 


Unfurling Acceleration 


An expression for the unfurling rate was given in the analysis section and can be 
written 


% = u = 



(Bl) 


where 


and 


or 


T = n sl K’ e 


sec 



e = 


X B " ^B 
L sl 


(B2) 

(0 < l B < L sl ) (B3) 

( L sl < l B < L p) <B4) 


An expression for the unfurling acceleration can be developed in terms of other known 
derivatives. This additional derivative was found to be necessary in order to minimize 
truncation error experienced with the use of a fifth-order Runge-Kutta numerical integra- 
tion technique. By assuming that the unfurling rate is a function of strain, the elastic 
modulus, and the linear mass density, the unfurling acceleration can be found by substi- 
tuting equation (B2) into equation (Bl) and differentiating: 
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APPENDIX B - Continued 


9u -L 9u / 8K sec dK sec 9u/8m' dm* 
u " 9e 9u/9e de a u /9e dZ B de J 


(B5) 


Evaluating the partial derivatives, equation (B5) can be expressed as 


n sl^sec i 

, , T ^sec 

1 (T + F re \ 

( Ib \ 

dm' 

2m'u 

n ,(K' de 

si v sec; 

m ' \ n slF-sec/ 

X B X B 
\ U " 

dZ B 


(B6) 


For practical systems, the term, 


T sec 

”sl( K sec) 2 d£ 

is of at least one order of magnitude less than 1. Therefore, equation (B6) can be 
reduced to 


• _ n sl K sec 
2m'u 


el 


1 /T + F re \ 
mf \ n sl K sec ) 




dm' 

dZ B 


(B7) 


where li is the unfurling acceleration, is the slope of the linear-mass-density 

d£ B 

curve, and e is the strain rate, which is given by the derivatives of equations (B3) 
and (B4); that is, 


or 


* = B*B “ ux b) 


e 


X B " u 
L sl 


(0 <lB< L s i) (B8) 

(l s i < £g < Lp) (B9) 


With the addition of equation (B7), the deployment motion is governed by seven first-order 
differential equations. 


Initial Conditions 

Initial values of the state variables of the vehicle and the deployment bag are well 
defined. However, since the initial values of the strain and the strain rate are undefined, 
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APPENDIX B - Continued 


the initial values of the unfurling velocity and acceleration are also undefined. Expres- 
sions can be developed for these values by applying the limit theory expressed by 
L'Hospital's rule. 

The strain during the initial unfurling sequence is defined as 


e = 



(BIO) 


Therefore, the limiting value of strain at the initiation of deployment can be expressed 


e 


o - 


AV 


- 1 


(Bll) 


where AV is the ejection velocity and the subscript o refers to the initial value. As 
shown previously, the suspension- line tension can be expressed in terms of the unfurling 
rate as 


T = m'u^ + F 


re 


(B12) 


or in terms of the strain as 
T = "slK' „* 


(B13) 


A cubic equation in the initial unfurling velocity can be found by equating equations (B12) 
and (B13), and evaluating the result at the initiation of deployment. The result can be 
written 

F. 


sec n . K’ 

u n 3 + — u n AV = 0 

m sl m sl 

Equation (B14) can be solved for the initial unfurling velocity. 

The strain rate during the initial unfurling sequence is given by 


(B14) 


a = A (Vb - ux b) 
6 b 


(B15) 

By successive application of L'Hospital's rule, the initial strain rate is determined to be 
u n x D ^ - u x D 

€ 0 = O B,o ^o B,o (B16) 

2u q 

where 

*B,o = ^v,o - v B>0 (BI7) 
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APPENDIX B — Continued 


The initial unfurling acceleration can be determined by substituting equation (B16) into 


equation 


(B7). /] 


Note that 


dm 


V 


dl 


B 


= 0. ] The result can be expressed as 


where 


and 


u o = x B,o 


4m',u 2 


(B18) 


AV t __sl_o 

u o K' 

sec 


AV = x 


B,o 


, ™ o 

min = 

sl n sl 


Discontinuities in Linear Mass Density 

Sketch B1 shows a typical profile of a precise linear mass distribution for a disk- 
gap-band parachute. The "spikes" in the linear mass density are due to cloth overlapping 



Sketch B1 
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APPENDIX B - Concluded 


and tape reinforcements in the canopy construction. These discontinuities in the linear 
mass density of the unfurling parachute result in a discontinuous unfurling rate. Recall 
the expression for the unfurling rate, which was expressed as 


u = 



where 


T n s iK sec e 


and 


or 


e 



B 



(B19) 

(B20) 

(0 < l B < L sl ) (B21) 


( L sl < l B < L p) ( B22 ) 


Since Xg and are continuous functions of time, the tension must be a continuous 
function of time. The unfurling resistance has been experimentally shown to have an 
essentially constant value during the unfurling of the suspension lines and another essen- 
tially constant value during canopy unfurling. Therefore, discontinuities in the linear 
mass density, along with this single discontinuity in the unfurling resistance, introduce 
discontinuities in the unfurling rate. As a result, in order to integrate the unfurling 
acceleration correctly, the integral must be reinitialized at each discontinuity. This can 
be easily done owing to the continuous nature of the tension. The value of the unfurling 
rate following a discontinuity in the linear mass density is given by 


u 2 



(B23) 


where subscripts 

1 denotes value of parameter prior to discontinuity 

2 denotes value of parameter following discontinuity 
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(a) Second balloon-launched flight test of Planetary Entry Parachute Program (ref. 6). 



(b) First flight test of Supersonic Planetary Entry Decelerator Program (ref. 7). 

Figure 3.- Comparison of calculated histories of unfurled length and flight test data for 

deployments of two disk-gap-band parachutes. 




(a) Second balloon launched flight test of Planetary Entry Parachute Program (ref. 6). 



(b) First flight test of Supersonic Planetary Entry Decelerator Program (ref. 7). 

Figure 4.- Comparison of calculated tension histories and flight test data for deployments 

of disk-gap-band parachutes. 
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(a) Precise linear mass distribution. 



(b) Averaged linear mass distribution. 

Figure 5.- Comparison of calculated deployment rate and unfurling rate, as influenced by 
different linear mass distributions, for the second balloon-launched flight test of the 
Planetary Entry Parachute Program (ref. 6). 
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(a) Precise linear mass distribution. 



(b) Averaged linear mass distribution. 

Figure 6.- Comparison of calculated deployment rate and unfurling rate, as influenced by 
different linear mass distributions, for the first flight test of the Supersonic Planetary 
Entry Decelerator Program (ref. 7). 
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